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r— I. ABSTRACT 

Context. Cool stars, companions to compact objects, are known to show Li abundances which are high compared to field stars of the 
Q , same spectral type, which are heavily Li depleted. This may be due either to Li production or Li preservation in theses systems. 

5_j . Aims. To measure the lithium isotopic ratio in the companion star of the neutron star X-ray binary Cen X-4. 

Methods. We use UVES spectra obtained in years 2000 and 2004 around the orbital quadratures. The spectra are analysed with 
' spectrum synthesis techniques and the errors estimated with Monte Carlo simulations. 

Results. We measure A(Li)=2.87 ± 0.20 and ^Li/'Li = 0. 12^q jjj at 68% confidence level. We also present updated system parameters 
with a refined determination of the orbital period and component masses i.e. 1.14 ± 0.45 Mq and 0.23 ± 0.10 Mq for the neutron star 
and companion, respectively. 

^ ' Conclusions. In our view the low level of *Li favours Li preservation scenarios, although Li production mechanisms cannot be ruled 

, out. In the case of preservation, no Li is freshly created in the binary, but the tidally-locked companion has preserved its original Li 

lO ■ by some mechanism, possibly inhibited destruction due to its fast rotation. 

m ■ 

^.f^ . Key words, accretion, accretion discs - binaries: close -stars: individual: Cen X-4 (=V822 Cen) - X-rays:binaries 

in 

J — 1. Introduction photometric modulations (Chevalier et al. 1989), and the deter- 

^ . , . ^ , . ,1, , . mination of the companion's radial velocity curve (Cowley et 

Spectroscopic analysis of the companion stars to X-ray bma- ^^gg. McChntock & Remillard 1990; Torres et al. 2002). 

_ > ries provides unique information on the nature of the accreting ^^j^^ ^-^^^^^^^ ^^^j^^j .^^^^^^^^ ^how a consistent velocity semi- 

_ compact objects (e.g.Casares 2001) and their formation history iitude(/r2 - 1 46 - 150 km s"') but a large scatter in systemic 

^ , ( sraehan et al 1999, Gonzalez Hernandez et al. 2004). The sub- velocities which tend to converge over the years from an inital 

^ . class of Soft X-ray Transients (SXTs) provide excellent testbeds 137.234 ^-1 (Rowley et al. 1988, McClintock & Remillard 

■ ■ ■ due to long periods of quiescence during which the companion ^^^^^ ^-i ^^^^^^^ 2002, D'Avanzo et al. 

star IS not overwhelmed by the (otherwise) mtense X-ray repro- 2005). The large systemic velocity and high galactic latitude 

cessed disc emission (van Paradijs & McC intock 1995). The ^ ^24°) is consistent with Cen X-4 being a member Galaxy 

study of SXTs has led to the discovery of a large population of ^^^^ population (Cowley et al. 1988). However, by integrat- 

black holes in the Galaxy, recognized by their large mass func- ■ ^. ^^^.^^ ^^^^^^-^ gravitational potential, 

tions (e g. Charles & Coe 2006, Casares 2005) exceeding 3 M^. Q^nzalez Hernandez et al. (2005a) propose that the binary was 

A sizable fraction of SXTs have shown Type-I X-ray bursts and, ^^^^^^ q^I^^^I^ 1^^^^ projected into its cuiTent posi- 

therefore they contain accreting neutron stars. ^.^^ ^ ^ ^^^^^ ^.^^ s^pemova explosion. The 

o^^^" ^ I ^ ^ °/ "^"T inclination angle has been constrained to / = 3 1 ° - 54° (Shahbaz 

SXp with two X-ray ouAursts observed in ^ et al. 1993) and the binary mass ratio to q = M2/MNS = 

,^0^^ ^u^^ (Kaluzienski, Holt & Swank 1980) q ^7 ^ q ^^^^^^ ^1 2002) which leads to Mns = 0.5 - 2.1 

V822 Cen, the optical counterpart of Cen X-4, was discovered ^ ^ q4 _ ^ ^^^^^^ ^^^^^.^^ 



and its companion. 



in the course of the 1979 outburst (Canizares, McClintock & 

Grindlay 1980). Since 1980 it has stayed in a quiescent level ^, . , , , , „^ 

. 10 o J u u »u » ^ e ■ ^ ■ The companion star has been classmed as a K3-K7 star 

at V^18.2 and it has been the target of intensive investigation , . , \ , . , , . , . 

f i« 1 u uv 1 ■ J »u u which must be substantially evolved in order to hll its ~ 

which led to the discovery of the 15.1 hr orbital period, through „^t^t^, . o. , 

; 0.6 Rq Roche lobe (see Shahbaz et al. 1993 and references 

Send offprint requests to: J. Casares therein). Detailed atmospheric parameters {T^s - 4500 + 100, 

* based on observations collected at ESO - Paranal in programmes \ogg - 3.9 + 0.3) and chemical abundances of the donor 

065-0447 and 073-0214 stai' were presented in Gonzalez Hernandez et al. (2005b). The 
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Table 1. Journal of Observations 



Date 


Wav. Range 


Number 


Exp. time 


Orbital phase 






of spectra 






25/4/2000 


4780-5755 

5837-6808 


10 


718 


0.09-0.22 


9/6/2000 


4780-5755 
5837-6808 


10 


718 


0.66-0.80 


3/4/2004 


3281^562 
4780-5755 
5837-6808 


4 


665 


0.64-0.68 


12/4/2004 


3281^562 
4780-5755 
5837-6808 


4 


665 


0.22-0.26 


16/4/2004 


3281-4562 
4780-5755 
5837-6808 


4 


665 


0.17-0.22 


21/4/2004 


3281-4562 
4780-5755 
5837-6808 


4 


665 


0.19-0.23 



low gravity provides additional support for an expanded atmo- 
sphere. Furthermore, a super solar metallicity [Fe/H]=0.23 + 
0.10 and moderately enhanced Ti and Ni is found, which sup- 
ports contamination from supernova products. The companion 
also presents an anomalous ^Li overabundance for a K type star 
(Martin et al. 1994; Gonzalez Hernandez et al. 2005b). This lat- 
ter peculiarity, which is also common to the companions of sev- 
eral black hole binaries, has been explained by spallation reac- 
tions in the course of the X-ray outbursts (Martin et al. 1994), 
although other mechanisms for Li production might be possible 
(e.g. Yi & Narayan 1997). Alternatively, the companion could 
have preserved its primordial Li abundance if the binary were 
younger than ^ 10^ yr or depletion efficiently inhibited by the 
action of fast rotation (see Maccarone, Jonker & Sils 2005). The 
origin of the large ^Li abundance could be addressed by measur- 
ing the isotopic ratio ^Li/^Li and, therefore, we have embarked 
in a project to obtain high resolution {R ~ 40000) spectroscopy 
of Cen X-4. 



2. Observations and reduction 

We obtained 20 spectra of Cen X-4 at the European Southern 
Observatory (ESO), OhservsAono Cerro Paranal, using the 8.2m 
Very Large Telescope (VLT), equipped with the UVES echelle 
spectrograph, on the nights of 25 April and 9 June 2000, cov- 
ering the spectral ranges /1/14780-5755 (hereafter green) and 
/il5837-6808 (hereafter red). A 1 arcsec slit was used resulting 
in a resolving power R=43,000. The exposure time was fixed to 
718 s and the observations were concentrated around the quadra- 
tures to minimize the effects of orbital smearing which, for the 
orbital parameters of Cen X-4, is in the range 0.2-10 km s~' i.e. 
smaller or comparable to our instrumental resolution of 7 km 
s-i. 

Sixteen more spectra of Cen X-4 were obtained on four 
nights in April 2004 also with UVES. The instrumental set- 
up was identical to the 2000 observations except that the 
DlCl dichroic was also employed, enabling us to obtain ad- 
ditional blue spectra covering /1/1328 1^562 (hereafter blue). 
Unfortunately, the count level was very low for some spectra and 
only 9 could be included in the analysis. The exposure time was 
set to 665 s which, for the phase of our observations, results in 
an orbital smearing < 7 km s~^ A full observing log is presented 
in Table 1. 



We extracted the 1 -dimensional spectra using optimal ex- 
traction within MIDAS ECHELLE context, and the pixel-to- 
wavelength caUbration was obtained with the observation of a 
Th-Ar lamp. We compared our reduction with the ESO pipeline 
and found that they are comparable, except for the 2000 data 
where we prefer ours because it provides better statistics, in par- 
ticular for the green spectra. All the spectra were flux caUbrated 
using nightly observations of ESO flux standards. We checked 
the accuracy of the wavelength calibration using the /16300.3 and 
/15577.4 [Ol] sky lines and found it to be precise within 0.4 km 
s"' for the two databases. 



3. Revised System Parameters 

Radial velocities of individual spectra were obtained through 
cross-correlating our red spectra (/1/15837-6808) with a synthetic 
template computed for the stellar parameters and metallicity de- 
rived by Gonzalez Hernandez et al (2005b), i.e. T^s = 4500, 
log^ = 3.9, [Fe/H]=0.23, using Kurucz LTE models (Kurucz 
1993). Prior to the cross-correlation, the synthetic spectrum was 
degraded to the resolution of our UVES data (7 km s"') and rota- 
tionally broadened by 44 km s ' to match the rotational velocity 
of the donor star in Cen X-4 (Torres et al. 2002). A spherical 
rotational profile with limb-darkening e = 0.65 was employed 
(Gray 1992). The main disc emission lines (H,i,, Hel /15875 and 
/16678) and interstellar/atmospheric absorptions (Nal D, /16300) 
were masked from the cross-correlation. A least-squares sine 
wave fit to our velocity points yields the foUowing orbital so- 
lution: 

7 = 189.6 ±0.2 kms"^ 
P = 0.6290522 ± 0.0000004 d 
To = 2451660.6432 + 0.0005 
K2 = 144.6 ± 0.3 kms'^ 

where Tq is defined as the inferior conjunction of the optical star. 
The quoted uncertainties are 68 percent confidence and we have 
rescaled the errors so that the minimum reduced is 1.0. Fig. 
1 presents our radial velocity points folded on the best fitted so- 
lution. We have also extracted velocities using our green spectra 
but the orbital elements obtained are identical except for a smaU 
shift of +2 km s~^ in the y-velocity. Note that the velocity semi- 
amplitude of the donor is extremely well constrained because 
we have concentrated our observations around the quadrature 
phases to minimize orbital smearing. Furthermore, the exquisite 
velocity accuracy of each data point (< 1 km s"^) coupled with 
our long baseline leads to a determination of the orbital period 
with unprecedented accuracy. 

We have also searched for evidence of period variations by 
comparing our Tq with the various spectroscopic zero phases re- 
ported in literature, which we list in Table 2. Given our long 
baseline we were able to estimate a zero phase for years 2000 
and 2004 independently. Fig. 2 shows the (O-C) diagram for 
To=2449163.934 (Torres et al. 2002) and assuming a constant 
period at P = 0.6290522 d. A quadratic fit provides a signifi- 
cantly better description of the (O-C) diagram than a linear fit 
(reduced of 4.9 vs 9.7) and suggests that the orbital period of 
Cen X-4 decreases in a timescale of Po/^* = (1.22+0.31)xl0''yr. 
Note, however, that this claim rehes on the very first epoch de- 
termination, obtained by Chevalier et al. (1989) through photo- 
metric lightcurves. If this is masked out a linear fit is marginally 
better and hence the evidence for a period variation vanishes. 
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Fig. 1. Radial velocities of CenX-4 folded on the orbital solution 
of Table 1 with best fitting sinusoid. Individual velocity errors 
are always smaller than 1 km s~' and are shown in the bottom 
panel. Open circles represent the 2000 data and solid circles the 
2004 data. 

Therefore, this result must be treated with caution and further 
observations are needed to confirm or otherwise a possible pe- 
riod variation in Cen X-4. Furthermore, a period decrease would 
be unexpected from evolutionary grounds: conservation of an- 
gular momentum leads to an increase in binary separation and 
mass transfer is expected to be sustained by nuclear evolution of 
the companion (King, Kolb & Burderi 1996). 

We subsequently produced spectral averages in the rest 
frame of the companion star for the 2000 and 2004 databases, 
using our best ephemeris. Following Marsh, Robinson & Wood 
(1994) we first computed the optimal Vrotsin; by subtracting 
broadened versions of our synthetic template (in steps of 1 km 
s"^) and minimizing the residual. We used a spherical rota- 
tional profile with linearized limb-darkening e = 0.65 and ob- 
tain Vrot sin; - 46^7 km s"' for both databases, slightly higher 
but consistent within Icr with Torres et al. (2002). The error 
is likely to be dominated by the uncertainty in e which is not 
known for the absorption lines (Collins & Truax 1995, Shahbaz 
2003). In an attempt to derive a more realistic error, we also 
computed the rotational broadening for the extreme cases where 
e=0-l and find Vrot sin j = 44-50 km s"^ or 47±3 km s'^ At 
this point we note that Vrot sin i is known to be modulated with 
orbital phase because of tidal distortion of the donor star (e.g. 
see Casares et al. 1996). It shows a double-humped oscilla- 
tion with peak-to-peak amplitude of ^10% and maxima at the 
quadratures, exactly when our data were taken. Therefore, our 
determination is likely to be overestimated by ^5% and hence 
we adopt Vrot sin / = 44 + 3 km s"' for the sake of the system 
parameters determination. This value of the rotational broaden- 
ing, combined with our refined K2, implies a binary mass ratio 
q = 0.20 ± 0.03 (Wade & Home 1988). We can now take our 
q value into the ellipsoidal solutions of Shahbaz et al. (1993) 
and get a much better constraint on the inclination angle of 
i = 33 - 45° (see their Fig. 3). This leads to a significant im- 



0.02 - 




-0.04 - 



-5000 5000 

N (cycles) 

Fig. 2. The (O-C) with respect to the linear ephemeris with 
P=0.6290522 d and To=2449163.934. Open circles mark zero 
phases obtained photometrically and soUd circles spectroscopi- 
caUy . The best quadratic fit is displayed, with the 1 cr uncer- 
tainty level marked in dashed lines. 



provementin the binary masses, i.e. M/^s = 1-14 ± 0.45 Mg and 
M2 = 0.23 + 0.10 Mq. The error budget is clearly dominated 
by the still large uncertainty in the inclination which needs to be 
refined through new high-quahty ellipsoidal fits. 

4. Spectral Distribution 

Optimal subtraction of our broadened template from the 2000 

and 2004 spectral averages yields relative contributions of the 
donor star to the total flux at /16400 of 78 and 89 percent, respec- 
tively. This implies that the accretion disc is dimmer in 2004, by 
a factor ^ 2.3, with respect to 2000. In addition, the mean EW 
of the Ha emission drops from ^ 40 A in 2000 to ^ 22 A in 
2004. This in turn implies that the emission line flux also de- 
creases by a factor ^2.1 between the two datasets. Significant 
variability in the accretion disc luminosity of Cen X-4 has been 
reported in previous works, both at short (i.e. minutes, Zurita, 
Casares & Shahbaz 2003) and secular (i.e. days. Chevalier et al. 
1989) timescales. Furthermore, Cen X-4 is one of the brightest 
SXTs in quiescence, with Lx ^ 4 x 10^^ erg s~' and important 
variability (Campana et al. 2004 and references therein). It is 
likely that the optical changes that we observe are triggered by 
reprocessing of the variable X-ray flux into the accretion disc. 

Figure 3 presents the average spectra of Cen X-4 in 2004, 
plotted in units and dereddened by E(B-V)=0.1 (Blair et al. 
1984). The three spectral ranges have been rebinned into 1000 
pixels and smoothed with a gaussian bandpass of 2 pixels. We 
also show our synthetic template, conveniently rescaled to 89 
percent of the total flux at /16400. The bottom panel displays 
the residual of the subtraction which represents the spectral dis- 
tribution of the accretion disc Ught. The continuum is well de- 
scribed by a standard oc A^" law with a = 1.4 + 0.1, in agree- 
ment with previous studies (Chevalier et al. 1988, Shahbaz et 
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Table 2. Observed Zero Phases 



Observed Xq 


NT 




Reference 


\ri.j LJ- z.'-t^yjyjyjyj ) 


(cycles) 






5778 3420(40) 


-5382 


-33.1 


Chevalier et al. (1989) 


6822.6121(90) 


-3722 


10.4 


Cowley etal. (1988) 


6909.4198(73) 


-3584 


8.9 


McClintock & Remillard (1990) 


8376.3664(90) 


-1252 


5.8 


Shahbaz et al. (1993) 


9163.9340(10) 








Torres et al. (2002) 


11660.6431(10) 


3969 


0.9 


VLT 2000 


12057.5761(10) 


4600 


2.0 


D'Avanzo et al. (2005) 


13099.9113(35) 


6257 


-2.3 


VLT 2004 





1 1 1 1 1 


1 1 1 1 1 


1 






Residual light 




H h- 



I I I l' I l' I I I 
4000 5000 6000 

WavclcngUi (A) 



Fig. 3. Top panel: Spectral distribution of Cen X-4, dereddened 
by Eb-y = 0.1, and synthetic template, scaled to 89 percent of 
the total flux at /16400. Bottom panel: spectral distribution of the 
accretion disc after subtracting the synthetic template from Cen 
X-4. 

al 1993). Furthermore, extrapolation to 2000 A yields an excel- 
lent agreement with HST UV fluxes reported by McClintock & 
Remillard (2000). The spectral distribution that we find is flat- 
ter than the standard a = 2.3 predicted for viscously heated 
optically thick accretion discs and suggests that irradiation of 
the outer disc regions may be important even in quiescence. 
On the other hand, the emission line fluxes yield the follow- 
ing Balmer decrement //„////j///-^///a=2.64/l/0.53/0.28, with 
F(Hii) = (2.14 + 0.06) X 10"'^ ergs cm^^ s^' A^'.This compares 
very well with the standard case B recombination for Tetj ^10"* 
K suggesting that, at the large densities typical of accretion discs 
(A^e ^ 10'^ cm"-^), emission lines are powered by photoion- 
ization (Drake & UMch 1980). The same result was also ob- 
tained for the quiescent black hole SXT V404 Cyg (Casares et 
al. 1993). This in turn supports a scenario where the optical vari- 
ability observed in both the continuum and line fluxes is trig- 
gered by reprocessing of the variable X-ray flux in the accretion 
disc. 

5. Li abundance and ^LiPU isotopic ratio 

5.1. Strategy of the analysis 

The D resonance doublets of ^Li and ^Li are separated by about 
0.16 A (the ''Li to the red), which implies that the ''Li Dl line 
is blended with the ^Li D2 line. In practice these form a unique 



unresolved feature even in warm slowly rotating stars, in which 
the feature is weak (see Smith et al. 1998 and Cayrel et al. 1999 
for examples). The companion star of Cen X-4 is rotating rapidly 
and, to further complicate matters, the line is strongly saturated. 

To analyse the Li doublet we proceed in a manner similar 
to that of Cayrel et al. (1999), who convincingly argued that 
the Li feature should be fitted with a synthetic spectrum using 
5 parameters: 1) placement of the continuum, 2) and 3) abun- 
dances of the two elements 4) macroscopic broadening of the 
lines and 5) wavelength zero point adjustment. The stars anal- 
ysed by Cayrel at al. were all slow rotators, thus the macroscopic 
broadening was dominated by instrumental profile and macro- 
scopic motions in the atmosphere, in our case it is dominated by 
the large rotation velocity. We thus used only Vrot sin / as fitting 
parameter, assuming a gaussian instrumental profile of 7 km s"' 
, as deduced from the Th lines in the calibration arc. Thus in our 
case the five parameters were: 1) placement of the continuum 
2) log(A?(''Li)/A?(Litot)) 3) A(Lif]4) y^t sin / 5) wavelength zero 
point adjustment. Of these 5 parameters 2 and 5 are correlated, 
quite obviously, however also 2 and 3 are strongly correlated, 
as we shall show below, due to the fact that the line is strongly 
saturated. 

5.2. Model atmosphere and synthetic spectra 

We used the same ATLAS LTE, plane parallel model atmosphere 
as in Gonzalez Hernandez et al. (2005b), with T^ff = 4500 K, 
logg -3.9 and [M/H] - +0.25. The Hthium doublet is heavily 
saturated, therefore we expect it to be seriously affected by mi- 
croturbulence. The determination of the microturbulent velocity 
requires the measurement of both weak lines (on the linear part 
of the curve of growth) and strong lines (on the flat part of the 
curve of growth). For a rapidly rotating star, like the companion 
to Cen X-4, the measurement of the weak lines is impossible, 
therefore one has no real handle on microturbulent velocity. The 
need for a microturbulent velocity arises from the limitations of 
our one-dimensional model-atmospheres and the need to take 
into account the velocity fields which characterize a real stellar 
atmosphere. If one were to use a full three dimensional model 
atmosphere there should be no need to introduce such a param- 
eter It has long be recognized that microturbulence is not a re- 
ally independent parameter, but can be calibrated as a function 
of effective temperature and surface gravity. One such calibra- 
tion has been derived for stars in the solar neighbourhood with 
similar metallicity and stellar parameters (Allende Prieto et al. 
2004), when applied to our adopted atmospheric parameters it 
provides ^ = 0.91 kms '. It may be questionable whether such 
a calibration is applicable to a stars which is rapidly rotating, 
since it has been derived for "normal" slowly rotating dwarfs 

' A(Li) = log(A'(Li)/A'(H)) + 12 
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Fig. 4. Best fit to the average spectrum of year 2000, after cor- 
recting every individual spectrum from the orbital solution pre- 
sented in Sect. 3. A veiling factor of 0.22 (defined as the relative 
contribution of the disc to the total flux) was corrected from the 
observed spectrum. 



Fig. 5. Comparison of different fits to the lithium doublet, with 
fixed zero shift. Solid line: our best fit with all other 4 parameters 
free which yields log(N(''Li)/N(Litot)) = -1.12. Dashed line: 
fit with log(N(''Li)/N(Litot)) = -3.5 fixed. Dashed-dotted line: 
same with log(N(''Li)/N(Litot)) - -0.47 fixed. The spectrum of 
CenX-4 was corrected from a veihng factor of 0.22. 



of spectral types K to F. One could however argue that merid- 
ional flows induced by rotation should be rather slow, for this 
rotational speed, and have a character of macrotrubulence rather 
than microturbulence, unless the flow becomes unstable, produc- 
ing small scale disturbances. Therefore we do not expect the mi- 
croturbulence to be too different from ~ 1 kms Our adopted 
value is therefore ^ - \ kms"', however we shall discuss the in- 
fluence of microturbulence on our results. The synthetic spectra 
were computed with the Linux version (Sbordone et al. 2005) 
of the SYNTHE suite (Kurucz 1993) We computed the emerg- 
ing specific intensity for a mesh of directions, then the emerging 
flux was obtained by integrating the specific intensity and taking 
into account the stellar rotation as described e.g. in Gray (1992) 
and implemented by the code ROTATE of Kurucz (1993). In 
this way limb-darkening is automatically taken into account. We 
computed a small grid of synthetic spectra with 3 isotopic ratios, 
3 Li abundances and 3 rotational velocities, namely A(Li)=2.35 
, 2.55 and 3.15, log(A?(''Li)/A^(Litot))=-4, -1.125 and -0.301 and 
Vrot sin /=40, 46 and 52 km s"' . We used the same line list as 
in Gonzalez Hernandez et al. (2005b) plus the complete list of 
Kurucz (1993) and the TiO bands from Schwenke (1998). For 
Li we took into account the full HFS structure listed by Kurucz 
(1995) as updated by Kurucz (2006). 

5.3. Line profile fitting 

We used as main data set the spectra observed in year 2000, 
since they have higher signal-to-noise ratio with respect to those 
observed in year 2004 (73 vs 34). Fits on the latter data pro- 
vide, nevertheless, consistent results. We unveiled our observed 
spectrum by 22 percent, as derived in Gonzalez Hernandez et al. 
(2005b) for the Li region. 

We wrote a fitting code, in which the minimum is sought 
numerically using MINUIT (James 1998) and the synthetic spec- 
tra were interpolated in our grid. One useful feature of MINUIT 
is that it allows to fix any of the parameters, if required. 

As a first try we fitted the whole region 6692.6-6723.5 A 
where the Li line is located, and our best fit to the coadded spec- 
trum of year 2000 is shown in Fig.|4] This fitting region is rather 
wide and we next attempted to fit only the region 6705.7-6709.4 
A and we find that the result does not depend on including 




Fig. 6. 3-D view of the histogram of Monte Carlo events with 
respect to Li isotopic ratio and total Li abundance for the model 
with ^ - I km s"' , Shift= km s ' (fixed) and veiling of 0.22. 



other lines except for the Li feature. In Table |3] we show the 
results obtained for different values of the veiling factor, the mi- 
crotrubulent velocity and leaving all 5 parameters free, or fixing 
at km s"' the velocity shift. 

The main motivation for searching for solution with a fixed 
kms"' velocity shift is that our coadded spectrum has been 
obtained by shifting at zero velocity all the individual spectra, 
to better than 0.1 kms"' , according to our error estimate, and 
hence we expect to find a zero velocity shift. To test the validity 
of a null velocity shift, we have performed the spectral fitting to 
the range /1/16300-6520 A and, indeed, obtain a zero point shift 
of kms"' . In the second place we want to test theories which 
predict production of ''Li in low mass X-ray binaries, a positive 
velocity shift will tend to conceal any possible presence of ''Li. 
Therefore a fit with a fixed km s"' velocity shift is useful and 
its result should be considered to be the "maximum" fraction of 
''Li compatible with the data. Finally also the results of Monte 
Carlo simulations (see next section) support a km s"' velocity 
shift as the most likely. 
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Table 3. Best fits for different initial parameters 



Veiling ^(kms-' ) Continuum Shift(kms-' ) K„, sin;(kms-' ) A(Li) log(A'(*Li)/yV(Lito,)) N(^U)/NCU) 



Single fits;;, 

0.22 1 1.07 0.36 49.6 2.97 -1.94 0.01 

0.22 1 1.07 49.6 2.87 -1.12 0.08 

0.28 1 1.07 49.3 3.12 -1.55 0.03 

0.22 2 1.07 50.1 2.78 -0.97 0.12 

0.22 2 1.07 49.3 3.05 -4 lO-"* 

Monte Carlo Simulations 

0.22 1 1.07 ±0.01 49 + 1 2.87 + 0.20 -0.97 + 0.19 0.12+°;;^ 

0.22 1 1.06 + 0.01 -0.42 + 0.62 49 + 1 2.83 +0.12 -1.05 + 0.22 O.lO+Sj^S 

0.22 2 1.07 + 0.01 50+ 1 2.77 + 0.17 -0.91 +0.22 O.M+jj^^ 

0.22 2 1.07 + 0.01 -0.38 + 0.50 50+ 1 2.73 +0.15 -0.88 + 0.48 0.15+55^1 

0.28 1 1.07 + 0.01 49 + 1 3.10 + 0.22 -1.01 +0.16 O-^^-om 

0.28 1 1.07 + 0.01 -0.77 + 0.70 49 + 1 3.01 +0.14 -1.02 + 0.22 0.10+!!i!5 



" Numbers in bold fonts correspond to parameters that have been fixed 



Inspection of Table |3]reveals that our derived Li abundance 
differs from that of Gonzalez Hernandez et al. (2005b). This is 
mainly due to the inclusion of the ''Li component but also, to 
some extent, to the readjustment of the continuum and rotational 
broadening, the different adopted microturbulence and spectrum 
synthesis code. In fact, when one fixes the Li isotopic ratio at -4 
in logarithmic scale (i.e. virtually no ^Li) the A(Li) obtained is 
almost identical to the value reported in Gonzalez Hernandez et 
al. (2005b). 

All the isotopic ratios provided in Table |3] imply 
very little ^Li (recall that the meteoritic isotopic ratio is 
log(A?(''Li)/A?(Litot)) = -1.123). The best fit, with 1 kms"' mi- 
croturbulence, 0.22 veiling and fixed kms"' velocity shift is 
equivalent to ^Li/^Li ^ 8% and is shown in Fig.|5]as a solid line. 
It is remarkable that in this case the ^Li content is consistent with 
the meteoritic abundance. 

In order to get a visual impression of the change the Li 
isotopic ratio makes to the profiles we also show in Fig. |5] 
two extreme fitting profiles: one with a very low ''Li content 
log(A^(''Li)/A^(Litot)) = -3.5 is shown as a dashed line, whereas 
another fit with very large ''Li content log(A^(^Li)/A^(Litot)) - 
-0.47 is plotted as dash-dotted line. These correspond to ^Li/'Li 
= 0.03% and 50% respectively. 

From Table [3] we may also deduce the impact on the derived 
Li abundance and isotopic ratio of the veiling factor and micro- 
turbulent velocity. We used the estimate of the veiling factor un- 
certainty by Gonzalez Hernandez et al. (2005b). For the reader's 
convenience the resulting uncertainties in A(Li), logarithmic ^Li 
fraction and ''Li/^Li ratio are provided in Table|4] which has been 
estimated for the case where the shift is fixed at km s"' . Since 
Gonzalez Hernandez et al. (2005b) assumed a microturbulent 
velocity of 2 kms ' , we have also quantified the effect of the 
microturbulence on the metallicity. We found that the metallic- 
ity decreases 0. 1 dex for an increase of 0.25 km s"' . This result 
is not surprising, since only strong Fe I lines could be analysed 
due to the large rotational velocity of the secondary star. 



Table 4. Impact of microturbulence and veiling factor on derived 
Li abundance and isotopic ratio 



Parameter 


AA(Li) ^\og(N(''U)/N(UiM)) 


A'('*Li)/yV(''Li) 


Single fits 


= + 1 km s ' 
Aveiling = +0.06 


-0.09 +0.15 
+0.25 -0.43 


+0.04 
-0.05 


Monte Carlo Simulations 


Af = -1- 1 km s ' 
Aveiling = +0.06 


-0.10 +0.06 
+0.23 -0.04 


+0.02 
-0.01 




2.2 2.4 2.6 2.8 3.0 3.2 

A(Li) 



Fig. 7. Histogram of Monte Carlo events with respect to total Li 
abundance for the model with ^ - I kms ' , Shift= kms ' 
and veiling of 0.22. 

5.4. Confidence intervals 

Although we have used fitting, it is clear that all the theorems 
which allow to derive errors and confidence limits from do 
not apply to the case of spectra, since adjacent pixels are always 



7 



0.0 




2.4 2.6 2.8 3.0 3.2 

A(Li) 

Fig. 8. Contour plot of the histogram of Monte Carlo events with 
respect to Li isotopic ratio and total Li abundance for the model 
with ^ = 1 kms"' , Shifts kms"' and veiling of 0.22. 



correlated, as pointed out by Cayrel et al. (1999) and Caffau et 
al. (2005). However, confidence limits can be reliably estimated 
using Monte Carlo techniques. To this end we performed Monte 
Carlo simulations by injecting Poisson noise, according to the 
standard deviations used in the fitting, into the best fitting syn- 
thetic spectrum and performing the fit with the same code used 
for the observed spectrum, for over 50000 samples. 

The results of four such simulations are given in the last six 
lines of Table [3] We would like to stress that the Monte Carlo 
simulations yield no significant dependence of the derived iso- 
topic ratios on the values adopted for the microturbulent velocity 
and veiling factor, and are all consistent within 1-cr. The two di- 
mensional histogram with respect to isotopic ratio and A(Li) of 
the simulation with ^ - I km s"' and a fixed km s"' velocity 
shift, is shown in Fig. |6l The strong correlation of the two pa- 
rameters is clear and results from the fact that the line is strongly 
saturated. The higher the ''Li fraction and the more the line is de- 
saturated, thus allowing to reach the observed equivalent width 
with a lower A(Li). In Fig. [8] a countour plot of the same his- 
togram is shown. From this distribution one may estimate the 
area which encloses roughly 68% of the Monte Carlo samples, 
which may be taken as akin to Icr for normal distributions. We 
find that the highest Li isotopic ratio compatible with the ob- 
servations at Icr is log(A^(''Li)/A^(Litot)) = -0.78 which corre- 
sponds to A^(<'Li)/A^(^Li) = 0.20. In fact 68% of the events are 
bounded by the following limits A?(''Li)/A^(Lito,) = 0.07 - 0.16 
and A(Li)=2.67-3.07. Therefore, we can claim a Icr detection of 
the Li isotopic ratio in Cen X-4 with N{^Li)/NCU) = 0.12+[; [;^ 
and A(Li) - 2.87 ± 0.20. The most likely value for the Li abun- 
dance and its uncertainty was infered from the histogram ex- 
tracted from the Monte Carlo simulations which is shown in Fig. 
|2l Note that, although the most probable value of A(Li) of the 
distribution coincides with the value found by the best single fit, 
the corresponding isotopic ratio is slightly different. 

It is interesting to note that in both the Monte Carlo simu- 
lations done leaving all five parameters free, the mean value of 
the velocity shift is consisten with kms"' , within 1 standard 
deviation. In our view this strongly argues in favour of fixing at 
Okms"' the velocity shift. 



6. Adopted A(Li) and isotopic ratio 

All the results of our fits and Monte Carlo simulations are listed 
in Table[3] whereas Table|4]summarises the impact of the veiling 
and microturbulent velocity. 

We have already discussed the issue of veiling factor and 
microturbulent velocity and from our point of view the best 
choice is 0.22 veiling factor and 1 kms"' microturbulent veloc- 
ity. Concerning the velocity shift we believe that the results of 
the Monte Carlo simulation and the process of optimal coaddi- 
tion of the spectra adopted, strongly suggest that this parameter 
should be fixed at km s"^ . We are left to chose between two 
possible couples of values for A(Li) and log(A^(*Li)/A^(Litot)): 
the best fit with a fixed kms"' velocity shift or the result 
of the corresponding Monte Carlo simulation. It is reassuring 
that the value of A(Li) is the same in both cases. The difference 
in isotopic ratio is small and boils down to ''Li/^Li = 8% for 
the direct fit or ^Li/^Li=12% from the Monte Carlo simulation. 
Our discussion would not change if we adopted one or the other. 
However we decided to adopt the result of the Monte Carlo sim- 
ulation. We believe that in a situation like this, with a minimum 
with a complex topology, and strongly correlated parameters, the 
Monte Carlo simulation does a better job in estimating the most 
likely parameters. 

We also note that < 4% of the events have A^(^Li)/A^(Litot) > 
0.4 which defines a 3-cr upper limit to the ^Li abundance of 
A^(''Li)/A^(^Li) < 0.66. 

7. Discussion 

The Li abundance in the companion of Cen X-4 is known to be 
high (Martin et al 1994, Gonzalez Hernandez et al. 2005b), and 
is confirmed by the present result of A(Li)= 2.87. The relevant 
aspect is that this value is higher of what expected for a cool 
K star, in which the convective zone is deep enough to destroy 
the surface Li. This characteristic is shared by other companions 
in X-ray binaries. Li detections have been reported in the com- 
panion stars to several SXTs, namely V404 Cyg (Martin et al. 
1992), A0620-00 (Marsh et al. 1994), GS 2000+25 (Filippenko, 
Matheson & Barth 1995) and Nova Mus (Martin et al. 1996). 

The surprising high abundance of Li observed in the cool 
companions of compact objects has been taken as an evidence 
of local Li production connected to the presence of relativis- 
tic particles around the compact objects. These are binary sys- 
tems in which a compact object, a neutron star or a black hole, 
is accreting matter from the companion and Li production can 
occur either in the accretion flow or on the surface of the com- 
panion, but through processes associated with the accretion (see 
Guessoum & Kazanas 1999 and references therein). Proposed 
physical mechanisms are a - a fusion in the accretion flow (Yi 
& Narayan 1997) or spallation of CNO elements by neutrons, 
on the star surface (Guessoum & Kazanas 1999). This leads to 
the prediction that other light nuclei such as *Li, Be and B must 
be synthesized besides ^Li, the amount depends on the details of 
the accretion process. 

However, Li in these systems is never found in excess of the 
presently cosmic, i.e. meteoritic value, A(Li)=3.3 which would 
be an unambiguous signature of fresh Li production. Thus it has 
also been suggested that instead of Li production we are deal- 
ing with the supression of the destruction mechanism, which is 
normally very efficient for late type stars. Theoretically rotation 
has been shown to counteract Li depletion and some evidence 
of reduced Li depletion has been found in tidally locked bina- 
ries in open clusters. Martin & Claret (1996) showed that high 
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Cen X-4 



Table 5. Li detection in SXTs 



System 


V f sin z 

' rot -^i*! » 


1 


' rot 


Spectral 


EWLi 


A (Li) 


Reference'' 




(kms"' ) 




(kms"^ ) 


Type 








T1 1 1 So-zlSn 

J 1 1 1 OT^'-rOU 


1 1 A -\-A 


f^9.-\-l 

DoitZ 




JV J - / V i 




<, 1 .o 




J0422+32 




45±2 


127±36 


M2-4V 


<0.21 


< 1.62 


(3,4) 


N. Oph 77 




70±10 


53+23 


K5-7IV 


<0.37 


< 2.96 


(5,6) 


GS2000+25 


86 ±8 


65±9 


95±11 


K3-6V 


0.150+0.090 


2.20 + 0.50 


(7,8) 


GUMus 


106+13 


54±2 


131 + 16 


K3-K4V 


0.420+0.060 


3.00 + 0.50 


(9,10) 


V616Mon 


93 ±4 


41±3 


142+11 


K2V 


0.245+0.030 


2.31 +0.21 


(11,12) 


V404 Cyg 


41 ±1 


55±4 


50+3 


KOIV 


0.290+0.030 


2.70 + 0.40 


(13,14) 


Cen X-4 


44 ±3 


40±7 


68+11 


K4IV 


0.440+0.025 


2.87'' + 0.20 


(15,16) 



" (1) Orosz 2001. (2) Gelino et al. 2006. (3) Harlaftis et al. 1999. (4) Gelino & Harrison 2003. (5) Harlaftis et al. 
1997. (6) Remillard et al. 1996. (7) Harlaftis et al. 1996. (8) Callanan et al. 1996. (9) Casares et al. 1997. (10) Gelino, 
Harrison & McNamara 2001 (11) Marsh, Robinson & Wood 1994. (12) Gelino, Harrison & Orosz 2001 (13) Casares 
& Charles 1994 (14) Shahbaz et al. 1994 (15) This paper (16) Shahbaz, Naylor & Charles 1993. 

* Note that in this case the Li abundance has been derived taking into account the presence of *Li, at variance with 
what done for the other stars. A fit to the data assuming no *Li would provide an abundance higher by 0.12 dex. 



enough rotational velocities may inhibit Li depletion in the at- 
mospheres of such late-type stars. On the contrary, Pinsonneault 
et al. (1990) have found that the effects of high rotational ve- 
locities on the Li depletion processes are small and these effects 
are not resposible for the Li abundance dispersion in Pleiades 
K-dwarfs with similar effective temperatures and masses, dis- 
carding a clear correlation between non-projected rotational ve- 
locities (obtained from photometric periods) and Li abundances 
(King et al. 2000). More recently, Maccarone, Jonker & Sills 
(2005) have argued that inhibited Li destruction, due to tidally 
locked rotation of the companion stars, could explain the large 
abundances observed in SXTs. X-ray binaries spend several Gyr 
of their lifetimes as tidally locked binaries and it is only when 
they come into the Roche-lobe contact that some Li depletion, 
through mass transfer, starts. Cataclysmic variables, which are 
systems related to the SXT, do not show the presence of Li. 
This fact has been taken as evidence of Li production in the SXT 
(Martm et al. 1995). However, Maccarone et al pointed out that 
cataclysmic variables have longer lifetimes compared with the 
systems ending with a black hole or a neutron star and that the 
companions in cataclysmic variables are not tidally locked for 
most of their evolution prior to the Roche lobe contact phase. 
These differences could explain the non detection of Li in the 
companions of cataclysmic variables. 

In Table 5 we list the EWs (corrected for veiling of the ac- 
cretion disc) and Li abundances of several SXTs, together with 
the stellar rotation Vmi, binary inclination (based on ellipsoidal 
model fits) and spectral type of their companion stars, as re- 
ported in the papers listed in column 8. A(Li) values are taken 
from the following works: Martm et al. (1996) for V404 Cyg, 
GU Mus and GS2000-H25; Gonzalez Hernandez et al. (2004) 
for A0620-00 and Gonzalez Hernandez et al. (2006) for XTE 
J 111 8+480. Upper limits to A(Li) have also been computed 
for GRO J0422-H32 and N. Oph 77 through the code MOOG 
(Sneden 1973) using the reported EWs and Teff of companion 
stars. Note that Li abundances in this table do not consider *Li, 
except for the case of Cen X-4. Viot is calculated from the mea- 
sured rotational broadening Viot sin / and the inclination. Table 
5 shows that there is no obvious correlation between A(Li) and 
Vrot but do seems to be a connection with the spectral type of the 
star: Li is only detected in companions with spectral type earlier 
than K5, suggesting that the depth of the convection layer seems 
to be an important factor This can be taken as circumstantial 
support to preservation models. 



A measurement to the ''Li/^Li ratio can, in principle, dis- 
criminate between production or preservation scenarios. This is 
because a + a fusion reactions occurring in the inner region 
of an accretion disc containing helium would produce signifi- 
cant amount of ^Li beside ^Li and ^Be which decays into ^Li 
with a lifetime of 76 days. Unfortunately, there are no compu- 
tations of a + a reactions in accretion discs with realistic (so- 
lar) abundances. Guessoum et al. (1997) examined the nuclear 
reactions that occur in the inner regions of accretion discs com- 
posed of pure helium. They computed the gamma-ray lines re- 
sulting from ^Li and ^Be deexcitation lines which produce a fea- 
ture at about 0.450 MeV and, incidentally, might have produced 
the 0.5 Mev emission feature observed in Nova Muscae. They 
conclude that the Li isotopic ratio is expected to be in the range 
0.1 < A^(^Li)/A^(^Li) < 10 for ion temperatures of the plasma 
varying from KT,=3MeV upto 50 MeV (Guessoum et al. 1997). 
Our analysis shows that A'(^Li)/A^(^Li) = 0.12+[5[5^ and IS con- 
sistent with the computations for ion temperatures ~ 4.3 MeV. 
Therefore, our value cannot rule out Li production models, how- 
ever we stress that the low level of ''Li is remarkable. This, to- 
gether with the possible correlation between A(Li) and spectral 
types and the non detection of a single X-ray binary with super- 
meteoritic A(Li), seems to support Li preservation scenarios. We 
also note that a new evolutionary picture has recently been pro- 
posed by Ivanova (2006) where LMXBs are powered by mass 
transfer from pre-main sequence donors. This would help to ex- 
plain the roughly primordial abundance of Li in Cen X-4 through 
youth and, hence, our value of the isotopic ratio. Clearly, in or- 
der to confirm our result, new computations of a + a fusion and 
spallation reactions in accretion discs with solar composition are 
required. 
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